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Random spin chains at quantum critical points exhibit an entanglement entropy between a seg- 
ment of length L and the rest of the chain that scales as log2 L with a universal coefficient. Since 
for pure quantum critical spin chains this coefficient is fixed by the central charge of the associated 
conformal field theory, the universal coefficient in the random case can be understood as an effective 
central charge. In this paper we calculate the entanglement entropy and effective central charge 
of the spin-f random Heisenberg model in its random- singlet phase and also at the critical point 
at which the Haldane phase breaks down. The latter is the first entanglement calculation for an 
infinite-randomness fixed point that is not in the random-singlet universality class. Our results are 
consistent with a c-theorem for flow between infinite-randomness fixed points. The formalism we 
use can be generally applied to calculation of quantities that depend on the RG history in s > f 
random Heisenberg chains. 
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I. INTRODUCTION 

Understanding universal behavior near quantum crit- 
ical points has been a major goal of condensed matter 
physics for at least thirty years. Quantum critical points 
describe continuous phase transitions at zero tempera- 
ture, where quantum-mechanical phase coherence exists 
even for the long-wavelength fluctuations that control the 
transition. Some quantum critical points can be under- 
stood via mapping to standard classical critical points 
in one higher dimension, but many of the most exper- 
imentally relevant quantum critical points do not seem 
to fall into this category. Furthermore, even quantum 
critical points that can be studied using the quantum-to- 
classical mapping have important universal features such 
as frequency-temperature scaling that do not appear at 
finitc-temperature critical points.^ 

One potentially universal feature of quantum critical 
points is the ground-state entanglement entropy, defined 
as the von Neumann entropy of the reduced density ma- 
trix created by a partition of the system into parts A and 
B: 



S = -Tr PA log2 PA = -Tr ps logj pb- 



(1) 



The entanglement entropy of the ground state at or near 
a quantum critical point can, in some cases, be under- 
stood via the quantum-to-classical mapping; an impor- 
tant example is that of quantum critical points in one 
dimension that become two-dimensional conformal field 
theories (CFTs), where the entanglement entropy in the 
quantum theory has a logarithmic divergence, whose co- 
efficient is connected to the central charge of the classical 
GET: 



2-5 



lim S^- log2 N, 



(2) 



where we consider A as a finite contiguous set of N 
spins and B is the complement of A in the chain. Away 



from criticality, the entanglement S is bounded above as 
N oo (the one-dimensional version of the "area law". ^) 
Surprisingly, the entanglement entropy, whose definition 
is closely tied to the lattice, is actually a universal prop- 
erty of the critical field theory, and hence independent of 
lattice details. At this time we lack a similarly complete 
understanding of critical points in higher dimensions; iso- 
lated solvable cases include free fermions^'*, higher di- 
mensional conformal field theories^'^, and one class of 
z = 2 quantum critical points.^" 

The connection between the central charge of GET's 
and their entanglement entropy implies that indeed for 
quantum critical points with classical analogs, the natu- 
ral measure of universal critical entropy in the quantum 
system (the entanglement entropy) is determined by the 
standard measure of critical entropy in the classical sys- 
tem (the central charge). In addition, it translates im- 
portant notions about the central charge to the realm 
of the universal quantum measure - the entanglement 
entropy. Zamolodchikov's c-theorem^ ^ states that the 
central charge c decreases along unitary renormalization- 
group (RG) flows. Therefore we conclude that the entan- 
glement entropy of GET's also decreases along RG flows. 
Stated this way, the strength of the c-theorem may apply 
to universal critical entropies in quantum systems that 
are not tractable by the quantum-to-classical mapping. 

One such class of systems is the strongly random ID 
chains with quantum critical points that can be studied 
by the real-space renormalization group (RSRG) tech- 
nique (see Ref. 12 for review). The archetype of such a 
system is the random Heisenberg s = 1/2 antiferromag- 
net, which exhibits the "random-singlet" (RS) phase. "'^^ 
In fact, disorder introduced to the pure Heisenberg model 
(which is a GET with c=l) is a relevant perturbation, 
which makes it flow to the RS phase. In Ref. 14 the 
disorder-averaged entanglement entropy is found to be 
logarithmically divergent as in the pure case, but with a 
different universal coefficient, which corresponds to an 
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effective central charge c = In 2, compared to c = 1 
for the pure case. This result has been verified numer- 
ically for the random-singlot phase of the XX and XXZ 
chains, which are expected to have the same critical prop- 
erties^^'i'^. 

The RS phase is but one example of an infinite ran- 
domness fixed point; special points which obey unique 
scaling laws. For instance, instead of energy-length scal- 
ing with 1/E ior pure quantum-critical scaling, we 
have:i3'i7 



Inl/S- L'^, 



(3) 



where in the RS phase, ip = 1/2. The infinite-randomness 
fixed points are, loosely speaking, the random analogs 
of pure CFT's, and therefore it is important to imder- 
stand all that we can about their special universal prop- 
erties, such as their entanglement entropy. Gencrically 
such points can be reached as instabilities to disorder 
of well-known CFT's (e.g. in the XX, Heisenberg, and 
transverse- field Ising model). Also, as we shall sec, ran- 
dom gapless systems exhibit RG flow between different 
infinite randomness fixed points. Another motivation for 
the study of the entanglement entropy in these systems 
is to understand if they exhibit any correspondence with 
the pure CFT c-theorem. Namely, does the entanglement 
entropy, or a related measure, also decreases along flow 
lines of systems with randomness? This question can 
be broken into two: (a) does the entanglement entropy 
decrease along flows between pure CFT's and infinite- 
randomness fixed points? (b) does the entanglement en- 
tropy decrease along flows between two different infinite 
randomness fixed points? 

The first of these questions was taken up by 
Santachiara,^* who showed that the random singlet phase 
of the parafermionic Potts model with n > 42 flavors 
has a higher entanglement entropy than the pure model. 
The second question, however, is still open. To answer 
it, and gain more insight into the cntangiomont entropy 
of random systems, we must consider systems which ex- 
hibit more complicated fixed points than the RS phase. 
In this paper we analyze entanglement entropy in the 
fixed points of the random Heisenberg spin-1 chain; par- 
ticularly at its critical point where the Haldane phase 
breaks down (the phase diagram of the spin-1 Heisen- 
berg chain is shown in Fig. 1), we refer to this critical 
point throughout as the Haldane-RS critical point. This 
extends our previous work by the authors and others 
on the entanglement entropy in a variety of "random- 
singlet" phases^'''^^"^° and opens the way for a similar 
calculation in the more complicated s > 1 Heisenberg 
chains. 

The strongly random critical points of s = 1,21-25 
s — 3/2,2^, and higher-spin^^ chains are ran- 
dom analogs of the integrable higher-spin Takhtajan- 
Babudjian chains. Consider s = 1 as an example: 
the Hamiltonian 



is critical, while without the bi-quadratic term it would 
be gapped; its critical theory'^" is referred to as the 
SU{2) k = 2 Wcss-Zumino-Witten model, with central 
charge c = 3/2. This critical theory has a relevant op- 
erator that corresponds to modifying the coefficient of 
the bi-quadratic term in the lattice model and thereby 
opening up a gap. There are similar unstable critical 
points in the phase diagram of higher-spin chains: al- 
though the generic higher-spin chain is cither gapped 
(for integer spin) or gapless (for half- integer spin), there 
is an integrable Hamiltonian, given by a polynomial in 
(Si • Si+i), that is gapless and critical with the central 
charge c — 3s/{s + 1). These critical points have addi- 
tional symmetry given by the SU{2) Kac-Moody algebra 
at level k = 2s. 

Before plunging into the details of our calculation, let 
us summarize our results. We find that the leading con- 
tribution to the entanglement entropy of the spin-1 ran- 
dom Heisenberg model at the Haldane-RS critical point 
is: 

5~ ice//log2L= — •(1.3327-10-3)-ln21og2L. (5) 

where the subtraction indicates the uncertainty in the 
results, which is an upper bound. The effective central- 
charge we find is thus: 



Jr.) 



1.232 



(6) 



Referring to question (a) above, c^y.^ < 3/2, i.e., is 
smaller than the central charge of the corresponding pure 
critical point. With respect to question (b), this effective 
central-charge is indeed bigger than the corresponding 
central-charge of both the Haldane phase, which van- 
ishes, and the spin-1 RS phase, which has: 

ln3 = 1.099. (7) 



„s=l RS 



H 



St 



Sj • Sj+i — (Sj • Si+i) 



(4) 



In our work we find that the methods used for the 
previously studied random-singlet-like critical points are 
inadequate to study the spin s > 1/2 more complicated 
critical points, where the history dependence of the renor- 
malization group is more complicated. The method de- 
veloped in this paper and applied to the spin-1 case pro- 
vides a general approach to the entanglement entropy 
of random critical points in one dimension accessible by 
real-space renormalization group. It also presents a well 
developed framework for other history dependent quan- 
tities, such as correlation functions and transport prop- 
erties. 

The following section reviews some basic results from 
the physics of infinite randomness fixed points, including 
the higher-symmetry points that exist in random spin 
chains. We review the strategy of our calculation in 
Sec. HI, and then proceed to derive the main results 
of this paper in sections IV- VII, where we obtain the 
renormalization-group history dependence for the s = 1 
chain and thereby its entanglement in a controlled ap- 
proximation and estimate the accuracy of the approxi- 
mation. The technique is compared to numerics for a re- 
lated simplified quantity of "reduced entanglement" for 
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which our technique provides exact results, with good 
agreement. In Sec. VIII we compare our results to the 
known value for a certain fine-tuned critical point of the 
pure spin-1 chain, which bears on question (a) above. 
We also compare the entanglement entropy at the two 
random fixed points of the spin-1 chain. 



II. REVIEW OF RANDOM-SINGLET AND 
HIGHER-SPIN INFINITE-RANDOMNESS FIXED 
POINTS 

A. Random-singlet phases - the simplest infinite 
randomness fixed points 

The spin- 1/2 Heisenberg chain is an antiferromagnetic 
chain at criticality, whose low-energy behavior is de- 
scribed by a conformal-ficld theory with central charge 
c = 1. Upon introduction of disorder, the low-energy be- 
havior of the chain flows to a different critical phase: the 
random-singlet phase. 

The random-singlet phase has very peculiar properties. 
The Hamiltonian of the system is: 



A: 



Ti — JjSi • Sj+i 



(8) 



Roughly speaking, the strongest bond in the chain, say, 

Ji, localizes a singlet between sites i and i + Quantum 
fluctuations induce a new term in the Hamiltonian which 
couples sites i — 1 and i + 2:^^'^^ 



Ji—\J-i 



i+1 



1,1-1-2 



2Ji 



Si-l • Si+2. 



(9) 



Eq. (9) is the Ma-Dasgupta rule for the renormalization 

of strong bonds. Repeating this process produces sin- 
glets at all length scales, as bonds renormalize into large 
objects. 

A useful parametrization of the couplings in the anal- 
ysis of the random spin-1/2 Heisenberg chain is; 



A = ln 



where Vt is the highest energy in the Hamiltonian: 
Q, = max,{ Jj}, 



(10) 



(11) 



and plays the role of a UV cnitofF. It is beneficial to also 
define a logarithmic RG flow parameter: 



r = in^ 



(12) 



where fio is an energy scale of the order of the maximum 
Jj in the bare Hamiltonian. In terms of these variables, 

and using the Ma-Dasgupta rule, Eq. (9), we can con- 
struct a flow equation for the distribution of couplings 



rfP(/3) _ dP{0) 
dv — a/3 



oc oc 



+P(0) / (i/3i / dM(p-p,-p,)P{l3i)Pm 





(13) 



Where the first term describes the reduction of fi, and 
the second term is the application the Ma-Dasgupta rule. 
For the sake of readability of equations, we denote the 
convolution with the cross sign: 



oo oo 







(14) 

Eq. (13) has a simple solution, found by Fisher, which 
is an attractor to essentially all initial conditions and 
distributions:^^ 



P(/3) = ^e-^^/^. 



(15) 



with X = 1- 

Many remarkable features of the random-singlet phase 
are direct results of the distribution in Eq. (15). In 
particular, its energy-length scaling, or, alternatively, the 
energy scale of a singlet with length L, is: 

L^~r = lnl/£; (16) 

with tp being a universal critical exponent: 

^ = 1/2. (17) 



B. Entanglement entropy in the random-singlet 

phase 

As mentioned above, in the random singlet phase of 
the spin-1/2 Heisenberg chain, the entanglement entropy 
of a segment of length L with the rest of the chain is:^* 



ln2 



■ log2 L. 



(18) 



The origin of this entanglement in the random singlet 
phase is simple to understand. Consider the borders of 
the segment L; the entanglement of Eq. (18) is due to 
singlets connecting the segment L to the rest of the chain. 
Each singlet contributes entanglement 1, and to calculate 
the total entanglement we need to count the number of 
singlets going over one of the borders of the segment L. 
Alternatively, the entanglement is twice the number of 
singlets shorter than L crossing a single partition (i.e., 
one of the barriers). 

In Ref. 14 we developed a method that allows to count 
this number of singlets as the RG progresses from high 
energy scales, in which short singlets form, up to the 
energy scales F ~ L^/^, where the singlets are of the 
same length as the segment length. Our method strongly 
resembles the calculation in Ref. 33. 
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In this paper wc will generalize this method in order to 
calculate the entanglement entropy of more complicated 
infinite-randomness fixed points in s > 1/2 spin-chains. 

The random-singlet phase can occur in the random 
Heiscnbcrg model of any spin s, but when s > 1/2, 
the randomness needs to be sufficiently strong, such that 
strong bonds restrict the two sites they connect to be in 
the complete singlet subspace. This implies that instead 
of entanglement 1, each singlet contributes: 

log2(2s + l) 

to the total entanglement. Therefore the entanglement 
of the spin-s random singlet phase is: 

ln2Jog2(2£+l) 
S r log2 L. (19) 

C. Higher spins infinite-randomness fixed points 

The random-singlet phase, with universal distribution 

(15) and length-energy scaling, Eq. (16), is but one ex- 
ample of an infinite randomness fixed point. In general, 
disordered systems may have a similar type of logarith- 
mic length-energy scaling relations, Eq. (16) with differ- 
ent ijj. This replaces the notion of the dynamical scaling 
exponent z in ~ 1/_E of pure critical points. Sim- 
ilarly, other infinite-randomness fixed points may exist 
with a different value of x in Eq- (15)) but with the same 
functional form of the distribution. The numbers tp and 
X are the critical exponents which parametrize different 
infinite-randomness universality classes. 

A prime example of different infinite-randomness criti- 
cal fixed points comes from Heisenberg models with spin 
s > 1/2. Following Refs. 21,22 and 26 which dealt with 
the spin-1 and spin-3/2 cases respectively, Damle and 
Huse have shown that a spin-s chain may exhibit infinite- 
randomness fixed points with: 

l/tp = 2s + l, X = 2s. (20) 

These fixed points were dubbed domain-wall symmetric 
fixed points for reasons that will become clear shortly. 



D. The spin-1 Heisenberg model 

This paper concentrates on the entanglement entropy 
of the random spin-1 Heisenberg chain at its critical 
point. The phase diagram of this system is shown in 
Fig. 1. Without disorder, the spin-1 chain has a gap 
(the "Haldane gap").^** The groimd state of the chain is 
adiabatically connected to the valence bond solid (VBS) 



of the AKLT chain, in which each site forms a spin-1 
singlet with both of its neighbors: (sec Fig. 2a) 

1^) = n - |0) (21) 

where and are the creation operators for the spin-up 
and spin-down Schwinger bosons. Each spin-1 site can be 
thought of consisting of two spin- 1/2 parts, symmetrized. 
In the VBS state it is those spin-1/2 parts that form 
the spin-1/2 singlet links over all bonds. These links are 
represented by the spin-1/2 singlet creation operator in 
terms of the Schwinger bosons: 

l=(arbl,-blal,). (22) 

As disorder increases, defects start occurring in the per- 
fect chain of singlet links (see Fig. 2b); these defects 
suppress the gap, until at a certain disorder ro the gap 
vanishes, and the chain enters into a Griffiths phase. ^^'^^ 
In this region there is still a line of singlets connecting side 
to side, but no gap. When a critical disorder Tc > rc is 
reached, the connection between the two sides also disap- 
pears, and the Haldane phase terminates. At disorders 
r > Tc the chain's ground state is the spin-1 random- 
singlet state. Our goal is to calculate the entanglement 
entropy at r = Tc- 

E. Entanglement entropy in the spin-1 Heisenberg 
model at the Haldane-RS critical point 

In order to calculate the entanglement entropy in the 
spin-1 Heisenberg model at its intermediate-randomness 
critical point, we employ the real-space RG technique. 
This technique generalizes simply to the spin-s > 1/2 
case. As in the spin-1/2 case, the first step is to find the 
strongest bond in the chain. But instead of putting the 
strongly interacting sites in a complete singlet, we just 
insert one Schwingcr-boson singlet (SBS) (Eq. 22) to the 
bonds' wave-function.^^ Effectively, this reduces the spin 
of each site by 1/2, and reduces the energy of the most 
excited state of the bond. If disorder is very weak, the 
SBS's form uniformly, and the AKLT state results. On 
the other hand, if the disorder is very strong, SBS's form 
in pairs between strongly interacting sites, and the spin-1 
random singlet phase results. 

Once more, the origin of the entanglement entropy is 
clear; each SBS going over a partition contribute entropy 
of the order 1 to the entanglement between the two sides 
of the partition. We need to count the entanglement of 
these SBS's until their length equals the size of the seg- 
ment under consideration. The challenge of this calcu- 
lation, however, is that each singlet going over the par- 
tition contributes an amount of entropy that depends 
on the singlet configuration forming both at higher and 
lower energies. This is due to the fact that each site 
can support two singlets connecting it to other sites, and 
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Haldane phase 



r RS Phase 

Lr C 



• • > > 

gapped * 

ungapped (Griffith phase) 



r 

(randomness) 



FIG. 1: Phase diagram of the spin-1 random Heisenberg modeL At r = the chain is in the gapped Haldane phase and 
its ground state resembles a valence-bond solid (VBS). As randomness is increased, the gap is destroyed at re, but the VBS 
structure survives up to the critical point, r = Tc- At r > the chain is in the spin-1 random- singlet phase. At the critical 
point, a different infinite randomness fixed point obtains, which has — 1/2, and x = 2. We concentrate on the entanglement 
entropy of this point. 




FIG. 2: (a) At very low disorder, the ground state of the spin-1 Heisenberg chain is well described as a valence-bond solid. 
Each spin-1 site is described by two spin-1/2 parts (black dots) that are symmetrized. Each site forms a spin-1/2 singlet to its 
right and to its left, (b) As disorder grows, defects appear in the VBS structure, and the gap is suppressed, (c) At very high 
disorder, a phase transition occurs to the spin-1 random singlet (RS) phase. 



therefore the ground state of the system is no longer a 
product state of site-pairs wave function. 

Another complication in the s = 1 case is the possi- 
bility of forming ferromagnetic bonds. This happens, for 
instance, if an even number of bonds localize one SBS 
each. The edges of this singlet chain have spin-1/2 each, 
which interact with each other ferromagnetically. At low 
energies in the RG process, this strong FM bond can be 
decimated, in which case the two spin-1/2 coalesce to a 
single spin-1. If the partition we are concerned with is be- 
tween these two spin-1/2, then after the FM decimation, 
the partition is lodged inside a spin-1 site. 

Thus to be able to calculate the entanglement entropy 
in s > 1/2 spin chains, we need not only a count of the 
SBS going over a partition, but a full knowledge of dis- 
connected singlet configurations: their probabilities, the 
energy scale at which they form, and their exact von- 
Neumann entanglement entropy. To carry out this cal- 
culation, we will use the domain-wall description^^'^^ of 
the Haldane- RS critical point. 



F. Domain- wall picture 

Our current understanding of the spin-1 critical point 
between the Haldane phase and the random-singlet phase 
relies on the Damle-Huse domain-wall picture. Let us 
first demonstrate this picture in terms of the spin-1/2 
random singlet phase. One can think of the random- 



singlet phase as forming through a competition between 
two possible singlet domains: Domain (1,0) with singlets 
appearing on odd bonds only, and domain (0,1) with sin- 
glets appearing on even bonds. The notation (a, 2s — a) 
with < a < 2s signifies a domain with a spin-1/2 singlet 
links on odd bonds, and 2s — a spin-1/2 singlet links on 
even bonds. This notation makes it easy to think about 
randomness as competition between different dimeriza- 
tions. For each domain, there is a probability pa to be 
of type (a, 2s — a), and also, for each domain, there is a 
transfer matrix, which tells the probability of domain a 
to be followed by domain a', which is Waa' ■ Note that: 

2s 

Waa' = 1. 

o'=0 

In the domain picture, at any finite temperature or en- 
ergy scale, the non- frozen degrees of freedom (i.e., spins 
that were not yet decimated) lie on domain walls. Thus 
in the domain wall between the (1,0) and (0,1) domains, 
there is one free spin-1/2 site (see Fig. 3a). This free spin 
interacts with similar spin-l/2's in neighboring domain 
walls through an interaction mediated by quantum fluc- 
tuations of the domain in between. Thus each domain 
of type a is associated with a bond between neighboring 
free spins, and has a distribution of coupling Pa{P), with 
f3 defined in Eq. (10). 

The renormalization of strong bonds is now described 
as the decimation of a domain. In the spin-1/2 chain. 
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whenever a domain is decimated its two neighboring do- 
mains, being identical, unite to form a single large do- 
main - a singlet appears over the domain, and connect 
the spins on the two domain walls. This is the Ma- 
Dasgupta decimation step. The random-singlet phase 
appears when the (1,0) and (0,1) domains have the same 
frequency. It is a critical point between the two possible 
dimerized phases associated with the two domains. 



These fractions are universal and a direct consequence of 
the bare spin of the model. 

The analysis of the domain-theory of the spin-1 critical 
point is explained in Ref. 27. The important aspects for 
the purpose of our calculations are the following: (1) All 
domains have the same bond strength distribution: 



±e-2/3/r 



(24) 



In s > 1/2 spin chains, the domain picture is richer. 
In the spin-1 Heisenberg chain there are three possible 
domains: (0,2), (1,1), and (2,0). The VBS is associated 
with the (1,1) phase, which has a uniform covering of the 
chain with spin- 1/2 singlet links. On the other hand, the 
strong randomness randoni-singlct phase in this system 
occurs when the competing domains are (2,0) and (0,2). 
This is completely analogous to the spin-1/2 random sin- 
glet phase, except that the domain walls consist of free 
spin-1 sites (see Fig. 3b). 

A general domain wall between domains a and a' can 
be easily shown to have an effective spin: 



(23) 



as each singlet link leaving the domain wall removes a 
spin-1/2 from it. 

Typically, the decimation of a domain involves form- 
ing as many singlet links as possible between the two do- 
main walls. If the two neighboring domain are identical, 
a' = a", then so arc the domain walls, and a full singlet 
is formed; this is the Ma-Dasgupta decimation rule in 
Eq. (9). If the two domain- walls are not-identical, and 
interact with each other anti-ferromagnetically, singlet 
links forming between the two domain walls will exhaust 
one of the doniain-wall spins, and the domain Da will be 
swallowed by the domain containing the exhausted spin. 

If the two domains neighboring a strong bond are dif- 
ferent, a' ^ a", and the interaction between the two 
domain-wall spins is ferromagnetic, the two spins unite 
into the domain wall between Da' and • For example, 
consider a (1,1) domain with an even number of links, it 
has to connect between a (2,0) domain and a (0,2) do- 
main; both domain walls will have spin-1/2. Upon dec- 
imation of the (1,1) domain, we are left with a domain 
wall between (2,0) and (0,2). which has a spin-1. 

Indeed at the critical point between the Haldane and 
random-singlet phases all three domains appear with 
equal probability - hence the designation - permutation 
symmetric critical point. Since each domain appears with 
the same frequency at the critical point, each possible do- 
main wall appears with the same frequency as well. Two 
domain walls: (0, 2) — (1, 1) and (2, 0) — (1, 1) are effec- 
tive spin- 1/2 's, whereas the third possible domain wall, 
(2,0) — (0,2) is a spin-1. Thus, at any finite but low 
temperature or energy scale, 2/3 of the unfrozen degrees 
of freedom are effectively spin-1/2, and 1/3 are spin-1. 



(2) each domain has an equal probability to be of any 
type: 



Pa 



1 



1 



2s 1 3 



(25) 



(this value indicates the probability after averaging over 
all possible neighboring domains), and (3) the transfer 
matrix Waa' is also the same for all a': 

Waa' = ^{1-Saa') = la-Saa')- (26) 
ZS 2 

Eqs. (24-26) give a complete description of the spin- 
1 VBS-to-RS critical point. From it we can deduce the 
energy-length scaling properties, Eq. (16): 



I.e., 



L^/^ ~ InOo/fi- 



(27) 



(28) 



III. 



OVERVIEW OF THE ENTROPY 
CALCULATION 



Combining the above results on this permutation sym- 
metric fixed point, in the next sections we calculate the 
entanglement entropy of the VBS-RS critical point. Un- 
like the spin-1/2 RS case, where each singlet contributes 
exactly entanglement 1, the spin-1 Heisenberg model can 
exhibit complicated singlet configuration consisting of 
overlapping singlet links, and even of ferromagnetic dec- 
imations. 

Nevertheless, the basic principle in our calculation re- 
mains similar to the spin-1/2 calculation outlined in Ref. 
14. We can quite generally write the entanglement en- 
tropy of a segment of length L as: 

1 In r 

~ -^Ceff l0g2 L = 2^—Stotal, (29) 

where, following Eq. (27), ~ i^/^ is the RG flow pa- 
rameter at the length scale L. Literally, this expression 
describes the accumulation of entropy due to configura- 
tions with average entropy StotaU which form on average 
when In F changes by This configurations connect the 
interior and of the segment to its exterior, on one of its 
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(0,1) 



1 ^" 



b. 



(1,1) 



(0,2) 



(2,0) 



FIG. 3: (a) Two domains are possible in the spin-1/2 Heisenberg model - (1,0) and (0,1). A domain wall between them gives 
rise to a spin 1/2 effective site, (b) In the case of a spin-1 chain, there are three possible domains: (1,1), (2,0), and (0,2), 
domain walls between them are effectively a spin-1/2 and spin-1 sites respectively. 



two sides, hence the factor of 2. The average entropy 
S total is: 



^ total 



(30) 



Prom Eq. (29) we can infer the effective central charge: 

'^Stotal 



Ceff = 



■ln2 



(31) 



Prom Eqs. (29) and (30) the challenge in the calcula- 
tion becomes clear. Wc need to find i, the probabilities 
Pc for each configuration, and the entanglement exhibited 
in each such configuration, Sc- 

In Sec. IV we develop a scheme that calculates the 
possible configurational histories in the RG process, and 
obtains £ as well as Pc for all configurations. In Sec.V we 
define and calculate exactly a 'reduced entanglement'. 
This is a simple and rough measure of entanglement that 
is determined just by singlet counting. In Sec. VP we 
calculate this measure numerically, thus verifying the an- 
alytical calculation in Sec. IV. In Sec. VI we calculate 
Sc of the various configurations. In Sec. VII we combine 
the pieces into the final answer. 



IV. CALCULATION OF THE ENTROPY 
ACCUMULATION HISTORY DEPENDENCE 

A. Approach to the history dependence 

As pointed out in the introduction, one of the compli- 
cations in the spin-1 entropy calculation is the fact that 
the entanglement of SBS's depends on the previous and 
also subsequent RG steps. This arises because SBS cre- 
ate correlations between partially decimated spins. Por 
this reason, to find the entanglement entropy we need to 
calculate the rate of formation of various finite-size cor- 
related structures. In the following we develop a system 
that follows the probability and energy scale of the these 
structures. 

Consider a partition across which we calculate the en- 
tanglement. Let us assume that we are at an interme- 
diate stage of the RG, in which domains are long, and 



correlations between domain-wall spins can be neglected 
(these correlations decay exponentially with the domain's 
length). Also, we assume that the partition-bond (the 
bond in which the partition is situated) was created by 
a Ma-Dasgupta decimation (see Eq. 9). As we shall see, 
the last assumption assures that entanglement from sin- 
glets that form in the ensuing RG steps is independent 
of the steps preceding the Ma-Dasgupta decimation. 

At this point, when the partition-bond is decimated, 
there are three scenarios: 

(1) If the partition-bond is anti- ferromagnetic (APM), 
and the two domains neighboring the bond are different, 
then one of these domains must be the (1,1) domain (oth- 
erwise we have a (2,0) and (0,2) domains surrounding a 
(1,1) domain which must be ferromagnetic as discussed 
below). A single SBS forms over the bond and the par- 
tition, and the bond is absorbed to the (1,1) domain. 
Only one domain-wall spin is decimated this way, and 
the other survives. The entanglement of the singlet just 
created depends on what happens next to the undeci- 
mated spin. 

(2) If the bond is ferromagnetic (FM), i.e., it corre- 
sponds to a (1,1) domain of even length between a (2,0) 
and (0,2) domains, the two spin-1/2 domain-wall spins 
coalesce to form a spin-1 effective site (a domain wall be- 
tween the (2,0) and (0,2) domains) which contains the 
partition. Needless to say the entanglement depends on 
what happens to the partition-site in later stages of the 
RG. 

(3) If the bond is AFM, but situated between two iden- 
tical domains, the two domain- wall spins connected by 
the bond are completely decimated. The entanglement 
entropy of this event can only depend on the previous 
decimation history of these domain walls. This is the case 
of the Ma-Dasgupta decimation. Indeed this decimation 
directly follows a Ma-Dasgupta decimation as assumed 
above, it contributes 5^ = 1 between two spin-l/2's or 
Se = log2 3 between two spin-1 's. 

Whereas the third scenario returns the partition-bond 
to the starting point of the discussion, the entropy of 
cases (1) and (2) above will be determined by the ensuing 
decimation of the sites near or at the partition until a Ma- 
Dasgupta decimation occurs. This is a general principle 
for all spin-S models; once the partition-bond is freshly 
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determined by a Ma-Dasgupta decimation, it means that 
all spin-excitations within the domain are gapped out, 
and quantum fluctuations above the gap give the sup- 
pressed Ma-Dasgupta coupHng. To count the entangle- 
ment entropy we need to follow the decimation history 
of the partition-bond starting right after a Ma-Dasgupta 
decimation, until the next one. 

Our calculation will follow the above lines; we will 
count how long it takes, in terms of the RG progres- 
sion, to form various configurations between two Ma- 
Dasgupta decimations. Since there are three possible do- 
mains in the spin-1 Heisenberg chain, we need to consider 
each possible domain separately. In Ref. 27, Damle and 
Huse show that at the spin-1 critical point, each of the 
three possible domains, (2,0), (0,2), (1,1), appear with 
the same probability. Thus, right after a Ma-Dasgupta 
decimation step, the bond over our partition, has proba- 
bility 1/3 of being each of the domains. Due to right-left 
reflection symmetry, the contributions due to the bond 
being domains (2,0) and (0,2) are the same. Therefore 
we only need to consider two possibilities. 

The main complication in our calculation is the pos- 
sibility of ferromagnetic decimation steps. A FM deci- 
mation renders our partition lodged inside a spin-1 site, 
which is also a domain wall. Our calculation is simplified 
by splitting the history analysis of the partition-bond to 
events from a Ma-Dasgupta decimation up to the for- 
mation of a FM partition-site, and events following the 
partition-site formation until a Ma-Dasgupta decimation. 
We start by analyzing the latter. 



C. Starting with a spin-1 site 

Let us now follow the decimation process from the 
point that the partition is lodged inside a spin-1. Wo de- 
note s{£) as the probability of the partition to be inside 
an untouched spin-1 site. Due to the additivity principal, 
Eq. (33), we can use the boundary condition: 



s{i = 0) = 1 



(34) 



and follow how this probability decays as the partition- 
site becomes correlated through SBS's to its environ- 
ment. 

The possible RG histories of a spin-1 cluster with dis- 
tinct entropy formation are shown in Fig. 4. Each of the 
eventualities in the figure has a probability distribution 
fn^ a function of £ — In P. Note that 
in the e„ and /„ eventualities the Ma-Dasgupta rule is 
applied and we get a new bond which contains the par- 
tition; at this point the spins we followed are singlet-ed 
out, and the entanglement of the configuration is unaf- 
fected by the subsequent RG stages. The s„ eventualities 
return us to the original state of this part of the calcula- 
tion, i.e., a partition contained within a site. The need 
to separate these eventualities by the subscript n, of the 
number of effective sites to the left (right) of the parti- 
tion is necessary since configurations with different n's 
produce different entanglement entropies. 

Given the scaling form of the bond-strength distribu- 
tion. 



B. Note on additivity and independence of RG 
times 

One of our goals is the calculation of the energy scales 
it takes uncorrelated configuration to form. As shown in 
Ref. 14, probability functions of configurations should be 
calculated as a function of 



^ = inr/ro. 



(32) 



In the next sections we will use the fact that this time 
£ is additive in the following sense: the RG time between 
the formation of a partition-site (i.e., a spin-1 site con- 
taining the partition) to a Ma-Dasgupta decimation into 
a partition-bond is independent of the RG history before 
the formation of the spin-1 partition-site. This is clear 
since all the information about the coupling strengths of 
the bonds that led to the formation of the partition-site 
is swallowed within the composite spin-1 site, and the en- 
suing decimations are only determined by bonds coupling 
the partition-site to the rest of the chain. Thus the RG- 
time it takes to go between Ma-Dasgupta decimations is 
the sum of the time for a partition-site to form, and for 
the partition-site to be decimated: 



t^MD-MD 



t^MD-PS 



£ps- 



MD 



(33) 



P(/3) = -exp(-2/3/r), 



(35) 



we can now calculate exactly the probability evolution of 
all eventualities. First, the probability Si will diminish 
due to decimations of bonds on either of its sides: 



ds{£) 
dT 



-2P(0)s(£). 



(36) 



and recall that £ and F are related by Eq. (32). 

Starting with e_i, the immediate-singlet eventuality: 



de-i{£) 
dT 



':SeP{0). 



(37) 



The interpretation of this equation is that the change in 
probability of the e_i eventuality, is P{0)dT - the proba- 
bility that the bond to the right (left) is decimated, times 
the probability that the domains to the left and right of 
the bond are the same, which is 1/2, times the proba- 
bility of the initial configuration, S£. Upon substituting 
P(0) we obtain: 



de-i{£ ) 
d£ 



(38) 



If the bond to the right (left) of the partition connects 

to a (1,1) domain, then the spin joins a long (1,1) domain 
to its right (left). This domain has coupling strength /? 
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^^+RtoL 




FIG. 4: The possible outcomes of a spin-1 site when our partition is located inside the site. When the bond to the right (left) 
of the site gets decimated, it either forms a complete singlet - in case the domain to the left of the bond is the same as that 
to its right - or makes the spin join a (1,1) domain, in which case the bond strength has probability distribution 5o(/3). These 
two possibilities occur with the same probability, 1/2, due to the domain-wall permutation symmetry. The (1,1) domain can 
be decimated by forming another singlet (eventuality eo), by forming a triplet (eventuality so) or by joining on its left (right) 
another (1,1) domain through a Ma-Dasgupta decimation (eventuality /o). It could also form another bond to its left, giving 
rise to the next line of diagrams, and the eventualities Si{l3), ei, si, /i. The subscript indicates the number of bonds to the 
left of the partition that are involved in the (1,1) domain. Note that the e„ and s„ eventualities occur on a partition-bond 
decimation (P — 0), when Dl = Dr and Dl 7^ Dr respectively, and thus with the same probability due to domain- wall 
permutation symmetry {Dr/i^ are the domains to right/left of the bond). 



between its walls, and hence this eventuality is charac- 
terized by the distribution: 6*0 (/3). This distribution can 
evolve under the RG due to decimations of other bonds to 
the right (left) of the partition, where the (1,1) domain is. 



But any decimation of its neighboring bonds will remove 
probabihty from 5*0 (/3), as will a Ma-Dasgupta decima- 
tion of the bond on which the partition sits, in the case 
of /3 = 0. Thus the flow equation for Sa{/3) is: 



J 



dSo{(3) ^ dSo{/3) 

dr dp 



-pmsom 



■m) X 5o(/32))-2P(0)5o(/3) 



-s,P(0)P(/3). 



(39) 



r 



The last term represents what happens when the bond 
to the right (left) of the partition site is decimated, and 
the second-nearest domain is (1,1), which happens with 
probability 1/2, hence the factor of 1/2. The second- 
to-last term describes the removal of probability due to 
decimations on either side of the partition-bond after it 
was formed. The brackets describe what happens when 
the bond to the right (left) of the partition-bond is deci- 
mated. The two factors of 1/2 in the brackets are the 
probability of the second-nearest domain on the right 
(left) being different than (1,1) in the first term, and 
(1,1) in the second term. The first term on the RHS de- 
scribes the reduction of the energy scale il. The cross 
denotes convolution, as defined in Eq. (14). 



But a decimation to the left (right) of the partition 
can have two outcomes - if the bond to the left (right) is 
connecting the (1,1) domain of the partition with another 
(1,1) domain to its left (right), then a decimation of this 
bond means applying the Ma-Dasgupta decimation, and 
we can stop following the RG; this is /o outcome, and its 
evolution is: 



dfoji) 

dr 



1 



P{0) / dpSoiP). 



(40) 



Alternatively, if the second-nearest domain to the left 
(right) is not (1,1), then a decimation of the bond to 
the left will create another spin-1/2 singlet over it, and 



10 



will lead to the distribution S'i(/?). Si{f3) has the same 
structure of the flow equation, except for the source term, 



J 



which is now fed from ^o, rather than s^: 



dS^ = + IPmsm + X SM) - 2P{0)Sm + ^P(0)5o(/3). 



(41) 



r 



again the factor of 1/2 in the source term is due to the 
probability of the second-nearest domain on the left not 

being (1,1). 

Going back to So{P), when the bond on the partition is 
strong, it gets decimated. If the two neighboring domains 
are the same, then we apply the Ma-Dasgupta rule once 
/? = 0, in which case we obtain the singlet-ed cluster cq. 
The evolution of this probability is: 



deoje) 
dT 



(42) 



where the factor of 1/2 is the probability that the two 
neighboring domains are the same. If they are different, 
a new spin-1 site is born, which is eventuality sq. It is 



easy to see that: 



dV 



dep 
dT ' 



(43) 



Note that we ignore what happens to the resulting new 
spin-1 configurations; subsequent RG of this configura- 
tions are the same as those considered in this section, 
and can be analyzed self-consistently. 

Eq. (43) completes the consideration of all first-level 
eventualities. It is easy to see that Eqs. (40) - (43) 
generalize to the case of starting with S'„(/3), which is the 
distribution of a bond with n singlets to its left (right), 
and a long (1,1) domain to its right (left). Thus the 
equations for /„, e„, s„ become: 



%a = ip(o)/d/35„(/3). 



ds^ = + iP(0)(5„(/3) +P(/3i) X S^m) - 2P(0)^„(/3) + iP(0)5„_i(/3), 



(44) 



dr 

dT 



^5„(0), 



r 



Before solving all equations, we note that the hardest 
part of solving Eqs. (40)-(44) is finding the distributions 
Sn{(3)- Here, an important simplification occurs: Eqs. 
(39) and (44) for Sn always admit a solution of the form: 

5„(/3) = a„P(/3). (45) 

Substituting Eq. (45) in the equation for Sn gives: 

^ = -4^ + iP(0)s,, 

OO 

Note that / rf/3S'„(/3) = a„. 



Starting with Eq. (36) it is easy to see that: 
ds£ dsi 



(46) 



de 



and: 



se 



-it 



(47) 



(48) 



e_i is immediately obtained as: 

1 



e-i 



(1 



(49) 



Note that this is the probability for a double singlet form- 
ing to the left, and the same probability applies to singlet- 
formation to the right. Hence the total probability of a 
double singlet forming is 1/2. 
Next is ao: 



e ao ~ J ^ 





At. 



(50) 



In fact, Eq. (46) is: 



e a„ = J die Un-i 





(51) 
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It is thus helpful to define: 



an 



(52) 



In terms of the a we have: 



ao = Jdil = e 



an = J dian-i = 



(53) 



From here we can find all other eventualities. The 
right-leaning singlet clusters e„ probabilities are: 



fn 



Jo 



die 



in+l 

(^^Tl)!' 



(54) 



for n > 0, which is extremely simple considering the tor- 
tuous way to here. The final bin probability to end up 
in, say, bin e„ is: 



f 

Jo 



die 



in+l 



{n + l)\ 4"+2" 
The sum of all e„ eventualities is: 



OO OO ^ 

^ e„ = ^ = — 



(55) 



(56) 



n=0 



n=0 



This is the probability to end up as a right-leaning sin- 
glet cluster. The same calculation applies to the s-bins 
and f-bins. Indeed - 6 x = 1/2, which completes the 
probability of 2e_i to 1. Note that the total probability 
of the spin-1 surviving and becoming a new spin-1 is the 
total Sn sum, which is also 1/6. 



The RG time of the process follows directly from the 
above discussion. Let us calculate the total time it takes 
for a spin-1 to be completely eliminated. This we do in a 
self consistent way: The time it takes to get completely 
decimated from the s„ outcomes, is the same as the total 
RG time we are seeking. Thus: 



4=2 /de_r£+2^ /"(d/„+de„)-£+2 ^ /"ds„-(£+4) 

(57) 

where the factor of 2 is due to the reflection symmetry 
of the configuration enumeration. Upon use of the defi- 
nitions of e„, /„, Sn we have: 



r /■OO /■O 

lh = 2 e-^^ide + 6 
6 Jo Jo 



-4e \ " 



en+2 



n=0 ^ ' 



^6 



1) = - + --- = - (58) 
8 9 16 12 ^ ' 



r 



Hence: 



(59) 



This is the RG time for a spin-1 partition-site to be dec- 
imated completely, and form an uncorrelated partition- 
bond. 



the bond being a (2,0) or (0,2) domain, or second, being 
a (1,1) domain. The latter is simpler, hence we start with 
the former. 

Right after a Ma-Dasgupta decimation, the bond dis- 
tribution R{(3) is quite different from P(/3). It is actually 
the convolution of P(/3) with itself: 



D. Getting from a (2,0) domain to a spin-1 site 

A more complicated analysis is reqiiircd for the cal- 
culation of the history of a partition-bond between a 
Ma-Dasgupta decimation and the formation of a spin-1 
partition-site. We must consider two possibilities: first. 



Qid) = P(/3i) X P(/32) 

= 7d/3i 7d/?2P(/3i)P(/32)%+/3,-/3) = ^/3e-2/3/r 



(60) 

As the RG progresses, this distribution may rebound to 
be closer to P(/3). The flow equation it obeys is: 



J 



« = « + 2.P,o,ip,A)JTO-P(oW). 



(61) 
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The second term represents the probabihty that a bond 
on either of the two sides (hence the factor 2) gets dec- 
imated (hence -P(O)). If the two bonds right next to 
the decimated bond are different (probabihty half), then 
nothing happens to the bond on the partition - the dec- 
imated bond either forms a new spin-1, or connects a 
(1,1) domain to a spin on the side of the partition. But 
if the domains near the decimated-side bonds are the 
same (probabihty 1/2 as wch), wc get a Ma-Dasgupta 
decimation that involves the partition bond, hence the 
convolution. In this case, we also need to remove the 
original probability from the distribution -R(/3), which is 
the origin of the last term. 



partition-bond be decimated, with the same domains on 
its sides. This invokes the Ma-Dasgupta rule, terminates 
this step, and begins a new cycle. The differential prob- 
ability of this happening is ^R{0)dT. But there are two 
sub-possibilities: the neighboring domains could be ei- 
ther (1,1), or (2,0)7(0,2), each with probability 1/2. The 
first feeds into p'_i, and the second to So we obtain: 



dp -I _ dp'_^ 

dr ^ dr 
12 2 i^-e 



4 r 3 



e " — e 



(66) 



Thus: 



The solution of Eq. (61) is straightforward, since R(f3) 
lives in the functional space spanned by P(/3) and Q{P). 
If we write: 



and: 



dp-i _ dp'_i _l ( -^i _ -4A 



(67) 



R{P)=a,P{f5)+hQ{P), 



we obtain: 



at 




(62) 



(63) 



The eigenvalues and eigenvectors of the matrix on the 
RHS are: 



-1 (1,1) 

-4 (2,-1) 



(64) 



Given that the initial conditions are a^=o = 0, 6^=o = 1) 
we obtain: 



ae = I (e-^ - e-«) 
be = l (2e-^ + e-«) 



(65) 



From here on, we can follow the same procedure as 
Sec. IV C. The first thing to consider is having the 



P-i 
and as ^ - 



5" 

oo, we have: 



■P-i 



P-i\e^oo =P-i\ 



1 
4 



(68) 



(69) 



So with probability 1/2, our bond gets decimated via a 
Ma-Dasgupta rule right away. 

The alternative possibilities arise from the case of the 
partition bond being decimated while its neighboring do- 
mains are different from each other - which happens with 
probability 1/2. In this case the partition bond joins a 
(1,1) domain to its right or to its left, each with another 
probability factor of 1/2. When this happens, we need to 
follow a distribution, So{P), which has the same mean- 
ing and flow equation as S'o(/3) from Sec. IV C, except 
for the source term, which is now jR{0)P{/3) (the P{/3) 
arises since it is the distribution of the (1,1) domain to 
which the partition-bond annexes): 



J 



dSoiP) dSo{f3) , 1 



dr 



dp 



-pmsoiP) 



P(/3i) X Som) - 2P{0)So{P) 



r 



-R{0)P{P) 



(70) 



As in Sec. IV C, once the partition enters the distribu- 
tion So (/3) three possibilities for termination of this stage 
exist: (1) form a singlet configuration via the e„ route, 
(2) be involved in a Ma-Dasgupta decimation with the 
bond to the opposite side of the (1,1) domain - /„ route, 
and (3) the partition-bond can be decimated while its 
neighboring domains are different, giving rise to a spin-1 
containing the partition - the s„ route. Also, S'o(/3) can 
flow into the additional 5'„(/3) just as before. 



It is easy to see that the flow equations for e„, /„, s„ 
and 5„(/3) are the same as Eqs. (44). The only difference 
is S'o(/3), due to its different source term. Let us solve 
So{P) for this case. Here too, we can write So{(3) = 
aa((')P(0)- Once we substitute this into Eq. (70), we 
obtain: 



da^ . 1 2 , 



(71) 
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R(|3) 




- ' p=o s„ (P) N - 

Sn 



FIG. 5: The possible outcomes of a partition bond of domain (2,0). When the partition-bond gets decimated, and it is between 
two identical domains, Dl = Dr, then the Ma-Dasgupta rule is applied. Depending on whether Dl ~ Dr — (1, 1) or not, 
eventualities p-i, p'_i are obtained. If Dl Dr, when the partition-bond gets decimated it joins a (1,1) domain, in which 
case the bond strength has probability distribution S'o(/3). The (1,1) domain can be decimated by forming another singlet 
(eventuality eo), by forming a triplet (eventuality sq) or by joining on its left (right) another (1,1) domain through a Ma- 
Dasgupta decimation (eventuality /o). It could also form another bond to its left, giving rise to the next line of diagrams, and 
the eventualities Si(/3), ei, si, /i. The subscript indicates the number of bonds to the left of the partition that are involved in 
the (1,1) domain. Note that the e„ and s„ eventualities occur when Dl ~ Dr and Dl 7^ Dr respectively, and thus with the 
same probability due to domain- wall permutation symmetry {Dr/l are the domains to right/left of the bond). 



As before, it is beneficial to follow uq = e^^a^. This 
obeys: 



(e--l). 



(72) 



The next step is to find a„, which are the probability 
'amplitudes' of Sn{P), as in Eq. (45). They obey Eq. 
(52). As it turns out, it is helpful to expand the expo- 
nent in Eq. (72) into a power-law. Upon integration we 
obtain: 



1 ^ 1 (3£)J+i 



from which we can immediately find: 



(73) 



- _ 1 1 (3^)J+"+^ _ I , u V (3^)' 
~ 3 ^ 3"+i (? +n+ 1)! ~ 3"+2 ^ j\ ' 

j=l j=0 

(74) 

From this result we can find e„ = /„ = s„ as a function 
of i: 



den 
'dl 



-u 



(75) 



and hence: 



en{(.) 



de 



1 



3"+2 



n+1 

E 



(My 
r 



\ J=o 

If we carry the integration to ^ — s- cx) we obtain: 

/ n+l -\ 

■3V 1 



^ 1 l l-(3/4)"+n 
\^ 4 1-3/4 J 



The sum over all e„ is just: 

1 1 



n=0 



1 



16 1-1/4 12 



(76) 



(77) 



(78) 



just as before. 

We are now ready to calculate RG times. The total RG 
time for this stage, between a Ma-Dasgupta decimation 
and the formation of a partition-site is: 
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As it turns out, we can carry out the sum over a„ directly: 



(79) 



And thus: 



oo 



n+1 



n=0 n=0 

C30 CX) 



lf>3« _ 1 
6^ 3 



E E 3^^ + 1/3 



j=0 n=j-l 

— Ip3i _ 1 _|_ 1 
6 2^ ^ 3' 



(80) 



4=2 1 / /I I 5 1 1^17 

° 3^ W' J \Q 2 37 8 3 8 12 





(81) 



r 



This is the average time it takes to become either a Ma- formation of a spin-1 partition-site (Fig. 6). 
Dasgupta decimated bond, or a spin-1 containing the par- 
tition (probabiUty 1/6). 



E. Getting from a (1,1) domain to a spin-1 site 

As in Sec. IV D, we first find the flow equation for 
A fresh partition-bond of a (1,1) domain can terminate the distribution i?(/3) of the partition-bond right after 
either through the formation of a full singlet, or by the decimation: 

I 



dR{P) _ dR{/3) 
dT ~ dp 



+ 2 • P(0) ( ^P(/3i) X R{p2) + Ir{0) 1 - 2P{O)Ri0). 



(82) 



r 



Note that there arc two differences with respect to the 
analog Eq. (61) in Sec. IV D. The brackets contain an- 
other term ii?(/3), which arises from the possibility of 
a bond to either side of the partition-bond being dec- 
imated, but with a different domain on its other side. 
This will extend the (1,1) domain of the partition-bond, 
but will maintain its strength. The other difference is 
the factor of 2 in the last term - this change is due to the 
fact that now any decimation of a bond neighboring the 
partition bond leads to its renormalization. In spite of 
these differences, the two additions cancel, and we obtain 
exactly the same result for R{l3) as in Eqs. (62) and (65). 

All that remains is to find the dependencies of the two 
probabilities e and s on £, where e is the probability 
of the partition-bond being decimated with its domains 
being the same on both sides, and ,s is the probability of 
being decimated with different domains surrounding the 



partition bond, and hence forming a spin-1 (see Fig. 6). 
It easy to see that these probabilities are equal, and that: 



de 1„,„, 12, 
— = -R(0) = --(e- 
dr 2 ^ ' T3^ 



(83) 



and hence: 



At 



)) (84) 



as ^ 00 we indeed obtain s(oo) = e(oo) = 1/2. 

The RG-time for this stage, i.e., between a Ma- 
Dasgupta decimation and the decimation of the partition 
bond (either FM or AFM), is: 



00 



dJl-i 



3 V 16/ 4 
(85) 
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FIG. 6: The possible outcomes of a partition bond of domain (1,1). When the partition-bond gets decimated, if it is between 
two identical domains, Dl ~ Da, then the Ma-Dasgupta rule is applied, otherwise, it forms a partition-spin-1. 



Note that now the probability of becoming a spin-1 (s 
eventuality) is 1/2. 



F. Total average time 

The average total RG time between two Ma-Dasgupta 
decimation is given as the sum of all the time of the 
limited histories times their probabilities. For instance, 
the total time it takes for a partition-bond which is a 
(1,1) domain to be completely eliminated is: 

1-513 

The last term is the product of the probability of forming 
a partition-site, and the time for this spin-1 to be com- 
pletely decimated. We use the results from Eqs. (59) and 
(85). 



(1,1) 



hi 



(86) 



A similar result applies to the (2,0) or (0,2) initial do- 
mains: 

^(.o)=V2)-^o + i4 = i^ + i4 = ^ (87) 

where we also use Eq. (81). Since this result is the same 
as for the (1,1) domain, the total average RG-time is: 



(■total — 



This number should be compared to the RG time be- 
tween Ma-Dasgupta decimations in the random-singlet 
phase, which is: 



nRS 
^total 



3. 



(89) 



Note also that this calculation is a generalized return-to- 
the-origin probability. 



V. REDUCED ENTROPY OF THE SPIN-1 
CHAIN 

After analyzing the history dependence and probabil- 
ities of the various singlet configurations, what remains 



is the analysis of the amount of quantum entanglement 
in them. Because of the complexity of the quantum en- 
tanglement of SBS configurations due to the correlations 
they produce, it is helpful to define a measure of entangle- 
ment which neglects these short range correlations. Such 
a measure would be a direct count of the number of sin- 
glets that connect the two parts of our chain, and can be 
calculated exactly. In this section we define this reduced 
entanglement entropy, E, analyze its properties, and cal- 
culate it both analytically, and in Sec. VF numerically, 
thus verifying our analytical calculation. 



A. Definition of the reduced entropy 

The ground state of the random spin-1 Heisenberg 
model is constructed through iterative decimations of 
strong AFM and FM bonds. The reduced entanglement 
with respect to a partition is defined as follows: An SBS 
connecting between one side of the partition to the other 
contributes reduced-entropy of 1. Some spin-1 sites are 
clusters that are the result of the decimation of a FM 
bond. If the partition is contained within such a cluster, 
a SBS that connects the partition site with the sites to 
its right contributes a fraction a to the reduced entropy, 
and an SBS to its left contributes a fraction 1 — a. 

a is supposed to capture the internal structure of the 
spin-1 partition-site, and in principle may vary from clus- 
ter to cluster. When two clusters combine, their weight a 
should be the average of their weights, cti2 = ^{'^i + '^2)- 
But using reflection symmetry of the problem, we now 
show that we can set 



1/2 



(90) 



without loss of generality. Every time that an SBS forms 
between a partition-spin-1 to a site on the right, we would 
have E = a per singlet. But because of the reflection 
symmetry of the problem we also need to consider the 
reflected configuration, which would have E — (1 — cr) 
per singlet. The average of these contributions is 1/2, 
independent of cr. Thus, for simplicity, we set a = 1/2. 
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Since the reduced entropy E neglects correlations be- The average contributions for this stage, together with 
tween the spin-l/2's degrees of freedom connected by the the partition-spin- 1 stage are: 



SBS's, it is a purely additive quantity: the reduced en- 
tanglement of a configuration is the sui 
entanglement of the individual singlets. 



tanglement of a configuration is the sum of the reduced En t\ = - ■ 1 + -Ei = - -\ — ^'^ (94) 

. - ■ ^^'^> 2 2 2 2 • 40 



B. Reduced entropy of a partition-site-1 



E. Average total reduced entropy, and reduced 
central charge 



Wc; start our analysis with a calculation of the average 
reduced entanglement of configurations that result from 
a partition site. Once we obtain our result for this inter- 
mediate stage wc will consider the average reduced en- 
tanglement of the beginning stages from a Ma-Dasgupta 
decimation up to the formation of the spin-1 partition 
site. 

The entropy contribution and probabilities of a 
partition-spin- 1 are illustrated in Fig. 7. The average 
contributions for this stage, and its futures are (going 
from right to left in the figure): 



Et = 

+ 1 . 1 + f 1 _ 1^ 



(91) 



where the last term self consistently adds the probability 
of forming a partition-spin- 1, which is 1/6, and it having 
the same entropy Ex. We get: 



Ex = 'l. 
^ 40 



(92) 



The average total reduced entropy is: 

Etotal = 5-^(1,1) + 1-^(2,0) 

— i -L _iZ_ -L II -L _iZ_ 

6 6-40 3-6 9-40 



(95) 



— 26 
18 



1 
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Wc are now in a position to calculate the reduced cen- 
tral charge of the spin-1 critical point. The entanglement 
of a segment L with the rest of the chain consists of: 



E = 2- ^ 



= 2-L^37|ia-lnLi/3 

14 (26 ]_\ 1 r 

3 3 V18 818/ ^^^-^ 



(96) 



Thus the reduced entropy central charge is: 



1 

12 



ln2 w 1.9171n2 



(97) 



which should be compared to Cr « 1.923 In 2 found nu- 
merically in the following section. 



Reduced entropy of a proton-bond of domain 
(2,0) 



F. Numerical evaluation of the reduced entropy 



The entropy contribution and probabilities of a 
partition-bond of a (2,0) domain are illustrated in Fig. 

8. The average contributions for this stage, together with 
the partition-spin- 1 stage are: 



E, 



2,0) 



i-l + i-2+i.2 + i.l + i.l + i£;i 



(93) 



= ^ + 

12 ^ 



47 

6-40 • 



D. 



Reduced entropy of a partition-bond of domain 
(1,1) 



The entropy contribution and probabilities of a 
partition-bond of domain(l,l) are illustrated in Fig. 9. 

Note that the original singlet of the (1,1) domain does 
not contribute, since we considered it when it formed. 



The previous sections have described a method for 
tracking history dependence in the RSRG in order to 
extract the universal part of the entanglement entropy. 
Since this method is fairly intricate, a direct numerical 
check on the results is worthwhile. Although the full en- 
tanglement entropy cannot be calculated in closed form 
but only in a controlled approximation, the "reduced en- 
tropy'' above can be calculated exactly. As a check, we 
now compute the reduced entropy by a numerical Monte 
Carlo simulation of the RSRG equations on a finite spin-1 
chain. 

For each subsystem size N, we average over 100 difi'er- 
ent intervals within each of 100 realizations. The initial 
distribution of Hcisenberg couplings is tuned to lie at 
the critical point of an infinite chain. Fig. 10 shows the 
result: the reduced entropy of an interval of size N is 
approximately 



S{N) = 1.9837 



1.923 In 2 



log^N. (98) 
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f—evenliialilies 
E=((1-<J)-KJ)=1 p=l/6 



s-cventikilUies 
E=O.5((l-o)+0)=O.5 



P=l/8 *: 



p=l/6-l/8=l/24 

E=(1-C!)+CT=1 



a=l/2 




Sj (P) 



p=l/6-l/8=l/24 
E=((l-o)+a)+l=2 



FIG. 7: Starting with a partition-spin- 1, the various possibilities are illustrated on the same chart as Fig. 4. a is the eflFective 
weight of the partition spin-1 that we begin with. Note that we omit here the future reduced entanglement in the case of the 
resulting partition-spin- 1. 

(2,0) 



f— eventualities 



p=l/6 





E=2 p=l/4 



E=l p=l/4 



i 



e— eventualities 




FIG. 8: Reduced entanglement starting with a (2,0) partition bond. The various possibilities are illustrated on the same chart 
as Fig. 5. Note that we omit here the future reduced entanglement in the case of the resulting partition-spin- 1. 



An overall statistical error in the numerical central charge 
of approximately 2% is expected. The agreement of 
the numerical value 1.923 In 2 with the analytic result 
(23/12) In 2 can be taken to verify the history depen- 
dence obtained above. It appears that systematic errors 
in the numerics resulting from finite interval size and fi- 
nite chain size (10* sites) are small. The remaining step is 
to use the same analysis of history dependence to obtain 
the true entanglement. 



VI. CALCULATION OF CONFIGURATIONAL 
ENTANGLEMENT IN THE SPIN-1 CRITICAL 
POINT 

The last step on our way to the entanglement entropy 
of the spin-1 chain at its critical point is to understand 
the quantum entanglement of each closed singlet configu- 
ration. Many such configurations arise in the decimation 
process, but fortunately, for a high accuracy evaluation of 
the entanglement entropy, only few classes of these con- 
figurations are necessary. In this section we will discuss 
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(1,1) 




E=0 p=l/2 E=l p=l/2 



FIG. 9: Reduced entanglement starting with a (1,1) partition bond. The various possibilities are illustrated on the same chart 
as Fig. 6. Note that we omit here the future reduced entanglement in the case of the resulting partition-spin- 1. 
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3.00 Lj-^'-' ^ ^ ^ ^ 1 — ^ — I 

10 100 

Interval size (sites) 



FIG. 10: Numerical calculation of the reduced entropy S'rod 
in the s = 1 chain for intervals of various sizes A*' within a 
chain of length 10^ sites. The fit is to 5'red = So + {c/3) logj N. 
The statistical standard deviation for each point shown is less 
than 5%. 



the entanglement in these configuration classes, both an- 
alytically and numerically. 



A. Simple configurations 

Let us first consider the configurations arising from a 
(2,0) domain, Fig. 5. The simplest eventualities involve a 
decimation of the (2,0) domain while its two neighboring 
domains are identical. There are two possibilities. In 
eventuality p_i the neighboring domains are (0,2), and 
the decimation of the partition-bond involves forming a 
complete singlet between the spin-1 sites to the left and 
right of the partition: 

\i'er)^l[aYbl~blaiy\0^r). (99) 



Since this is a Ma-Dasgupta decimation, the entangle- 
ment due to this step is independent of future steps. The 
entanglement between the left and the right sites is sim- 
ply: 

5p_, =log2 3 (100) 

since a trace of the left (or right) sites lead to an SO (3) 
symmetric density matrix for the right (left) site, and the 
spin-1 is completely undetermined, hence all three states 
contribute equally to the entropy. 

The second possibility is that the (2,0) partition-bond 
is near (1,1) domains, which is a p'_i eventuality. In this 
case, the domain walls surrounding the partition-bond 
are spin-1/2, and the decimation gives rise to one SBS, 
which gives: 

V, = 1- (101) 

A spin-1/2 degree of freedom is determined by the SBS 
on either side of the partition. 

If the partition bond gets decimated while it is sur- 
rounded by two different domains, then the bond joins 
the (1,1) domain by the formation of an SBS. Since the 
(1,1) domain is assumed to be very large in the scaling 
limit, the domain walls on both its sides are uncorrelated. 
Nevertheless, we know that the partition is at a small fi- 
nite distance to one of the (1,1) domain edges. In the 
next stages of the RG, the (1,1) domain expands, until 
it gets decimated in one of three ways. / eventualities 
will result in the (1,1) domain of the bond joining an- 
other (1,1) domain through a Ma-Dasgupta decimation 
on its short side. This results in a closed cluster, with 
the entanglement entropy: 

Sf = 1. (102) 

A more difficult configuration occurs if the (1,1) do- 
main containing the partition bond is decimated between 
two identical neighboring domains, while the partition 
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n 



c. 





FIG. 11: (a) A typical e-eventuality configuration, (b) Using 
the chain-rule singlet row reduction, we can reduce a chain of 
n + 1 spin-1 sites with two dangling spin-l/2's, to just two 
exact spin-l/2's, which are related to the dangling edge spin- 
1/2's as in Eq. (105). (c) Application of the chain-rule to the 
configuration in (a). 



bond is a distance n from the edge of the (1,1) domain - 
e„ eventuality of Fig. 5. In the absence of any ferromag- 
netic decimations, the entanglement of such a cluster can 
be calculated exactly, using the S0{3) symmetry, and an 
important simplification. This simplification serves as a 
key ingredient of the numerical calculation of the entan- 
glement entropy of more complicated configurations, that 
do involve ferromagnetic decimations. This is described 
in the next section. 



B. Singlet- row reduction 

To simplify the understanding of configurations such 
as in Fig. 11a, we will derive a rule that will replace 
the Schwinger-boson singlet chains which are completely 
contained in one side of the partition. Consider a chain 
of n -I- 1 spin-1 sites, whose state consists of a singlet- 
row, i.e., a single connects sites m and m + 1 for all 
TO = 0, . . . ,n— 1. Sites and n have a dangling spin-1/2, 
thus there are four singlet-row states, which we write as: 



IT (1, 1) T> = "n Ubl+i - bl^l+i) |0> 

IT (1,1) i) = 4^rn (al^I+i - ^l^I+i) io> 

li (1,1) ^)=blaiu(^rbli~M+i)\0) 

li (1, 1) i) = hXh\ "n - io) . 



(103) 



The dangling spin-l/2's at the edge are eventually locked 
into singlets, or ferromagnetic bonds which stretch past 
the partition. The sites i = 1, ...,n— 1 in the singlet 
row are inert, and do not contribute to entanglement, 
therefore calculation of the entanglement with them in- 
cluded is going to be quite difficult since they increase 
the Hilbert space involved by a factor of 3"""'^. They do 
have an important role, however, which is to facilitate 
the correlations along the singlet chain. 

But the correlations of the singlet-row can be taken 
into account while removing the interior inert sites. The 
correlations carried in the states in Eq. (103) are exhib- 
ited by the fact that the four states are not orthogonal. 
In particular: 



(T (1,1) iii (1,1) T) = (-ir 



(104) 



This overlap is easily found, since it is the product of the 
amplitude of the state ml = for all sites < i < n. 
Note that the states in Eq. (103) are not normalized: 

(T (1,1) TIT (1,1) T) = (i (1,1) ill (1,1) i)=iVn 
(T (1,1) ilT (1,1) i) = (i (1,1) Tli (1,1) T)=iVu 

(105) 

Let us now define four new states that will be mutually 
orthogonal: |TT) , |T-i),UT), Iii)- We can capture the 
correlations in the singlet-row by writing: 



IT (1,1) T) = v^lTT> 
li (1,1) i) = VW^lii) 
li (1,1) T) =«liT) + MTi) 
IT (1,1) i) -a|Ti) + feUT> 

where a and h are real numbers which obey: 



(106) 



+ 9 ^ iVi 
2ab= {-ly 



Ti, 



(107) 



(108) 



which is solved by: 

a = i (v/A^Ti + (-l)" + - (-1)") , 

b = k {^N^ + (-1)" - ^N^ - (-1)") . 
The norms Ni^i are found in App. A to be: 

A'n.lp- + (-!)■■) ■ 

JVtl = i (3"" - (-1)") ■ 

The eigenvalues of the reduced density matrix of states 
that consist of singlet-rows written just using the dan- 
gling spins, as in Eq. (106), will be the same as those 
for the reduced density matrix involving all inert spin-1 
sites, since the inner product of the states in Eq. (103) 
and (106) are the same. This presents an extreme simpli- 
fication in terms of numerical evaluation of these eigen- 
values. 
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C. Entanglement of the level-0 (no FM clusters) 
e-eventualities and f-eventualities 



Using the reduction of singlet-rows as in Eq. (106) we 
can reduce the configuration in Fig. 11a to that of Fig. 
11c. The entanglement entropy of the e„ eventuality is 
then easily calculated. 

The first step is to write the state in Fig. 11a while 
ignoring the correlations along the singlet-row: 

|^) = (|Tl)|i4)-Ul)|T4))(|T2)U3)-U2)|T3)) 

= |Tl)|i4)|T2)|i3)-|Tl)U4)|i2)|T3) (110) 
-Ul)|T4)|T2)U3) + Ul)|T4)|i2)|T3). 

Sites 1 and 2 are actually representatives of the dangling 
spin-l/2's in the singlct-row, and therefore we apply to 
them the assignment of Eq. (106): 

1^)^ VWltl>U4)|T2>U3) 
-a|Tl>U4)U2)|T3)-foUl>IT4>U2>|T3) 
-aUl>|T4)|T2)U3)-MTl)U4)|T2)|i3) 

+y^Ui)IT4)U2)|T3) 

with a and b defined as in Eq. (108). The entanglement 



(111) 



entropy follows from the eigenvalues of the reduced den- 
sity matrix of the state j-^) over sites 3 and 4. From 
SO (3) symmetry we can argue that the form of the re- 
duced density matrix is: 

Pi2 = tr34p = 

(ITT) (TTI + lii) (iil + I ITi + iT) (Ti + iTI) 
+/3ilTi-iT)(Ti-iTI 

(112) 

where this representation breaks the reduced density ma- 
trix into the triplet and singlet sectors respectively, where 
it is already diagonal. We also have a = 1 — /?. To find a 
and (3 all we need is to know the matrix clement of one 
of the triplet states in pi2. This is easily found: 



1 



7V- 



TT 



1 



TT 



2N^^ + 2a2 + 2P 2 iV|| -|- N^i 



(113) 



where the denominator is due to the normalization of the 
wave function in Eq. (Ill), and we used Eq. (107). This 
is then: 



a 



i-iy 



(114) 



4 4-3" 

The entanglement entropy of this configuration is then: 



= I - (, - ,„Mi - ») ^ - (2. t£) (i . i^) - (1 - tlL) (1 _ <_!£) 



The entanglement for the f-eventualities is remarkably 

simple, since in this case there is a spin-1/2 singlet cross- 
ing the partition, but this singlet is uncorrclated with the 
domain walls on the left and right of the partition bond. 
The configurational entanglement of an f-eventuality is: 

4? = 1- (116) 



D. Entropy of level-1 eventualities neglecting 
history dependence of cluster 

An approximate analytical answer for the entangle- 
ment can be obtained using the self-consistency method 
already at level-1, i.e., considering configurations with at 
most one FM decimations, by assuming that whenever a 
spin-1 cluster forms, its structure corresponds to of 
Figs. 4 and 5. This implies that the two spin-l/2's par- 
ticipating in the ferromagnetic cluster are uncorrelated 
with the partition bond, and are infinitely removed from 
it. In this case the entropy contribution from level-n and 



level-n-|-l become independent, and the entropy can be 

summed self-consistently, as described in Sec. VII A. 

For this purpose we need to derive the configurational 
entanglement of e„ and f eventualities of spin-1 sites, as 
in Fig. 4, assuming the initial FM spin-1 cluster consists 
of two uncorrclated spin-1/2 sites. Let us begin with 
the simpler f-evcntuality. This case is depicted in Fig. 
12a. The entanglement of this configuration, with the 
spin-1 partition-site containing two spin-l/2's which are 
uncorrelated before they form the FM cluster, is easily 
computed to be: 

5« = 2-ilog23. (117) 

Next we find the entanglement of the e„ eventualities of 
this setup, depicted in Fig. 12b. This setup is simplified 
to a six-site configuration, with the partition situated in 
the middle (Fig. 12c). Once the three spins to the right of 
the partition are traced out, we are left with a 3-spin-l/2 
density matrix, which for brevity we do not quote here, 
of a complicated mixed state. Due to rotational SO (3) 
symmetry, however, this density matrix breaks down to 
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a. 



b. 



c. 




FIG. 12: (a) A typical f-eventuality following a FM decima- 
tion. The entanglement entropy across the partition-site can 
be calculated exactly if we neglect correlations between the 
two spin-1/2 components of the partition site, a state which 
is denoted by Soo (b) A typical e-eventuality configuration af- 
ter a FM decimation, (c) Using the row-reduction rule we 
reduce the configuration from a to one of only six spin-l/2's. 
Neglecting internal correlations in the partition site, we can 
calculate this entanglement entropy exactly. 



J 



I3n = 



In = 



The above formulas, Eqs. (120) and (121) obtain for 
n > 1. The cases of n — —1, n — 0, and n = 1 need 
to be calculated separately, but their contribution to the 
entanglement entropy is very simple. The cases of e_i 
and eo implies the spin-1 partition site is decimated via 
two singlets that link it to two uncorrelated spin- 1/2 's 
to one of its sides. Upon tracing out of that side of the 
partition, we obtain a reduced matrix that describes a 
spin-1/2 with equal probability of pointing up or down, 



invariant sectors corresponding to the decomposition: 

1/2x1/2x1/2 = 3/2 + 1/2 + 1/2. (118) 

One sector is the spin-3/2 subspace, and two other in- 
variant subspaces will correspond to the total spin of 
the three sites being 1/2. Furthermore the restric- 
tion of the reduced density matrix on each of the in- 
variant subspaces, should be of diagonal form. In 
the basis |3/2,m)(TO = 3/2, . . . , -3/2), |l/2, m) (to = 
1/2,-1/2), |1/2',to) (to = 1/2,-1/2) the reduced den- 
sity matrix is: 

yOi,2,3 = diag{Q; In: In} (119) 

Finding a, /3, and 7 can be done analytically, a is the 
easiest to obtain, since we know from rotational symme- 
try that the state |tiT2T3) is an eigenstate. We readily 
obtain: 

«n = ^(l-(-l/3)") (120) 

Obtaining /3„ is more involved, but can be 
done by looking at the restriction of the re- 
duced density matrix to the invariant subspace of 
= 1/2 (U1T2T2), |Tii2T3), |TiT2i3)), and re- 
moving from it the subspace of |3/2,TO2 = 1/2) = 
75 (U1T2T2) + |Tii2T3) + |TiT2i3)). This leaves us with 
an easily diagnosable 2x2 matrix with eigenvalues: 



(121) 



hence, eg^'' = e-^\ = 1. 

An analysis of the ei case can proceed along similar 
lines to the one above and to the analysis of the e„ even- 
tualities as in Eq. (114), to obtain: 

4'^=4.4log2l8 + 2.-^log2^ (122) 



^ (^5 + (-1/3)" + y^(4 + 2 • (-1/3)")^ - 24 • (-l/3)« + 24/3^ 
+ (-1/3)" - y^(4 + 2 • (-1/3)")' - 24 • (-1/3)" + 24/32 



Thus the entanglement of cLii cluster contSiining r spin-l pa-rtition-site isi 

1 n = -l,l 

5W = <{ 4.^1og2l8 + 2.^1og2f n=l (123) 

-4a„ log2 - 2/3„ log2 /?„ - 27„ logj 7„ n > 1. 
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VII. ENTANGLEMENT ENTROPY OF THE 
SPIN-1 CHAIN 

Wc arc now at a position to combine the above results 
into the entanglement entropy, as given in Eqs. (29) and 
(30), to obtain: 



Ceff = — = — ln2. 



In Sec. IV we found that: 



e = 



(124) 



(125) 



Let us first find the average entanglement with the 
partition bond being a (2,0) domain initially. Looking at 
Fig. 5 we have: 

where Sc is the average entanglement entropy due to 
configuration c with c — e, f, s. The entanglement due 
to the level-0 e-eventualities starting with a (2,0) or (0,2) 
domain is given by: 



-q(O) 
•^(2,0) e 



00 ^ 

^ PpU ■ 1 +PpU ■ log2 3 + Yl ■ (130) 



n=0 



In the remains of this section, we will extract the value 
of 



S total — 'y^^PcSc 



(126) 



where Pc is the probability for a configuration c to occur, 
and Sc is the entanglement entropy associated with it. 

One major obstacle is presented by the ferromagnetic 
decimations, which give rise to complex quantum states. 
A decimated configuration can involve any number, n, of 
spin-1 ferromagnetic decimation steps, but with a rapidly 
decaying probability of 1/6" (see Eq. 56). The configu- 
rational entropy Stotai can be obtained to arbitrary accu- 
racy by calculating exactly the configurational-entropy of 
all configurations with n -~ 1 ferromagnetic decimations, 
and applying the approximation that the n-th level Fer- 
romagnetic decimation is between two spin-l/2's that are 
infinitely far from the partition, and thus correspond to 
the Soo eventuality. This provides a result for Stotai which 
has accuracy of: 



AS 



total 



6-6" 



(127) 



In what follows we obtain the computer- free n = 1 result, 
and using Mathematica to calculate and sum the config- 
urational entropy, we go up to level n = 3, with accuracy 
of ~ 10-3. 



where Pp'_^ = Pp_^ = 1/4. 

The entanglement due to level-0 f-entanglement start- 
ing with a (2,0) domain is: 



^(0) 

'(2,0) / 



= E4 

n=0 



(0) 



8-4" 



(131) 



To complete the calculation we need to add to these the 
entanglement of the spin-1 cluster that forms due to the s- 

eventualities, 5'(2|o) s- this section we approximate the 
spin-1 cluster as consisting of uncorrelated spins. This 
produces average configurational entropy we denote as 
5g , and thus: 



(132) 



where the one in the brackets is due to the SBS formed 
on the partition-bond in the process. 

The equivalent contributions to levcl-0 when the 
partition-bond is a (1,1) domain initially is quite sim- 
ple. With probability 1/2 the partition bond undergoes 
a Ma-Dasgupta decimation which has configurational en- 
tropy 1, which implies: 



0(0) 



(i:i)e = 1/2- 1 = 1/2. 



(133) 



A. Entanglement entropy from level and 1 



The remaining probability is that the partition-bond is 
decimated ferromagnetically. Thus: 



The entropy of the configurations computed exactly in 
Sec. VI C and VI D can be summed up to give the n = 1 
approximation for Stotai- In this approximation we have: 



Stotai — 2'S'(2,0) + 



(128) 



where 5'(2.o) and S'(i,i) are the average entanglement en- 
tropies across partition-bonds that are initially (2,0) and 
(1,1) domains respectively. 



= S(i 



(1,1) e 



-^(0) 
"^(1,1 



(1,1) s 



(134) 



We now calculate S^ . Looking at Fig. 4, we see that: 

(135) 



sT = si'^ + sf^ + si'^. 



The entanglement of an e„ cluster containing a spin-1 
partition-site is: 
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-^(1) _ 5 1 
-^6 00-3-^+32^0 



(1) 



^8-4 



- (-4a„ log2 an - 2/3„ log2 /3„ - 27„ log2 7„) . 



(136) 



where we make use of the configurational entanglement 
wc find in Sec. VI D, and where the Q!„, /3„, 7„ are de- 
fined in Eqs. (120) and (121). 

iSj^'' is much simpler, and from Eq. (117), which gives 
the configurational entropy, we get: 



-(1) 

'/ oo 



logs 3). 



(137) 



After considering all e- and f- eventualities, only one 
more possibility remains, which is the formation of a 
new spin-1 cluster via a ferromagnetic decimation. If we 
once more pursue the approximation of assuming that the 
spin-1 partition-site forms out of two spin-l/2's that are 
infinitely far from the partition, then wc can self consis- 
tently connect the entanglement of this eventuality with 
the one calculated in this section. Thus we write: 



-o(l) 



^(1) 



6 ' 



(138) 



Isolating , and plugging in all numbers from Eqs. 
(117) - (136) we obtain the approximate cluster entropy 
due to a spin-1 partition site: 



1.33495. 



(139) 



Note that this is an overestimate since we neglect corre- 
lations within the spin-1 effective site. These correlations 
reduce the entanglement entropy. 

Combining the results above, and substituting into Eq. 
(128) we obtain that the average configurational entan- 
glement is: 



Stotai = ^-1.16748+^ 



1.48271 = 1.3776. 



(140) 



where the tt^ describe what happens to the cluster at 
each level of FM decimation. A scries can terminate at 
■Ko = P-i, p'_i, en , fn , or begin by forming a FM 

spin-1 partition site with ttq = Sn • The following p-1 
steps can also be any sequence of FM spin-1 site for- 

R/L 



mation, tt^ 

_ R/L 
— ^n^ 



The last step, TTp, can be any of 
, fnj^, which renders the cluster decimated. 



Note that at levels p> 0, e„ has rip 



-1, 0, 1,.... 



Hence, an exact calculation can be written as: 



Stotai — 



E E 

P=0{7ri|i=0, 



^'^i 1 *^('7ro07ri07r20...07rp 



i=0 



(142) 

We know the probability of each eventuality exactly, as 
explained in Sec. (IV). The remaining part is the config- 
urational entropy. As explained above, we can calculate 
this entropy analytically for simple configurations. But 
to obtain a conclusive answer, we carry out a higher level 
calculation using the computer program Mathematica. 
Using the reduction trick of Sec. VI B, we can simplify 
any configuration of level-n to a density matrix involving 
at most 4 -F 3n spins. This allows an exact numerical 
computation. 

We carried out this computation up to level-3 config- 
urations. The result for the average configurational en- 
tropy was: 



5tota! = 1.3327 -0.001. 



(143) 



Note that this answer is an overestimate (hence the sub- 
traction sign for error indication), and thus can not be 
taken to be exactly 4/3. The closeness of the answer to 
4/3, however, is quite mysterious. 



B. Entanglement entropy from level 3 self 
consistent approximation 

In principle, we could obtain an exact result for the 
configurational entanglement by summing up the entan- 
glement of each configuration times its probabilities. The 
probabilities of each and every configuration arc well 
known, but at levels higher than one FM decimation, 
we can only obtain the configurational entropy using a 
computer. Nevertheless, this part of the calculation can 
be automated. 

The exact solution can be expressed using the following 
pattern. We can encode the eventualities of decimated 
configurations as: 



TTo O TTi O 7r2 O . . . O TTp 



(141) 



C. Effective central charge results 

The effective central charge of the critical spin-1 chain 
is thus: 



^eff 



-• 1.3327- ln2 = 1.232 
3 



(144) 



and the entanglement entropy between a segment of 
length L and the rest of the chain is: 



- 1.232 logaL. 

where we use the level-3 result above. Note that c 



(145) 

e// 



16/9 — e where e w 10~^. In the next section we discuss 
the significance of our results in the context of infinite- 
randomness fixed points. 
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VIII. DISCUSSION 

Several random quantum critical points in ID are now 

known to have logarithmic divergences of entanglement 
entropy with universal coefficients, as in the pure case. So 
far, all analyzed systems were infinite randomness critical 
points in the Random-singlet universality class. ^^'^^'^^'^^ 
Summarizing the results for these systems is possible with 
the formula: 

5~ ^InDlogsL, (146) 

where D is the dimension of the gapless sector of the 
Hilbert-space on each site, which generalizes to the quan- 
tum dimension in the case of the golden-chain. -"^^ Hence 
the 'effective central charge' for these critical points is 
simply: 

Ceff = In D. (147) 

The present work investigates the disorder-averaged 
entanglement entropy of a different universality class 
within the infinite-randomness framework: the spin-1 
Haldane-RS critical point. We indeed find a different 
result than Eq. (147), which arises from a more complex 
structure of the spin-1 Haldane-RS point. We find the 
entanglement entropy for this fixed point to be: 

S^^Ceff log2 ^ = y • 1-3327 • In 2 log2 L. (148) 

Our calculation required both the exact probabilistic de- 
scription of the low-energy structure given by RSRG, and 
the individual determination of the quantum entangle- 
ment of each spin configuration in the low-energy struc- 
ture. 

An open question is whether other physical properties 
are related to the universal coefhcient. In pure ID quan- 
tum critical systems described by 2D classical theories 
with conformal invariance, the central charge c controls 
several important physical properties beyond the entan- 
glement entropy. An example is the specific heat at low 
temperatures: because the central charge determines the 
change in free energy when the 2D classical system is 
compactified in one direction (i.e., put on a cylinder), 
the quantum system at finite temperature has a contri- 
bution to the specific heat that is proportional to c. 

The central charge also gives some information about 
the renormalization-group flows connecting different crit- 
ical points because of Zamolodchikov's c-theorem:^^ RG 
flows are from high central charge to low central charge, 
compatible with the heuristic view of the RG as "inte- 
grating out" degrees of freedom. Ref. 38 goes as far as 
defining a c-function for a finite-temperature quantum 
system, which decreases with decreasing temperatures, 
as it flows to a stable conformally-invariant fixed point. 
It is logical to ask"^'^** whether entanglement entropy can 
similarly be used to predict the direction of real-space 
RG flows. As pointed out in the introduction, Sec. I, 
this question separates in two: 



1 . flows from pure conformally-invariant points to infinite 

randomness points, 

2. fiow between two different infinite-randomness fixed 
points. 

Recently, a counterexample for question (1) was found 

by Santachiara:"'^* for the random quantum parafermionic 
Potts model with iV > 41, the random critical point was 
found to have larger entanglement than the pure critical 
point, even though randomness is relevant at the pure 
critical point. This point too obeys Eq. (147) with D 
N. Until now, there were no models in which we could 
ask the second question; the spin-1 random Heisenberg 
model provides the first test of case (2). 

The critical point of the spin-1 Heisenberg model we 
found the effective central charge 

16 

c'-J/j = 1.232 w y In 2. (149) 

As shown in Fig. 1, this point is unstable towards the 
spin-1 RS point with entanglement entropy: 

cf/y. = In 3 = 1.099 < c'-J'/j . (150) 

On the weak randomness side, the critical point is unsta- 
ble towards the Haldane phase (with topological order, 
but with a suppressed gap due to Griffiths effects. ^^■'^^ 
The Haldane phase does not have a InL term at all. 
Hence, within the critical points of the random Heisen- 
berg model, effective c does decrease along fiows, which 
agrees with case 2 above. 

The Haldane-RS critical point of the random spin-1 
chain is thought to be the terminus of a flow begin- 
ning with the pure spin-1 chain in Eq. (4), which is a 
k = 2 WZW theory with central charge c = 3/2. Since 

c),jI- — 1.232 < 3/2, the effective Central charge indeed 
also decreases along the flow line from the pure to the 
random critical fixed point, as posited in case 1 above. 

In order to make clear the similarity between the 
random-singlet phase of the spin-1 chain and the random- 
singlet phase of the spin- 1/2 chain, it seems worthwhile to 
introduce a modified central charge that is the coefficient 
of the 1/3 log2 L divided by the measure of the local un- 
gapped Hilbert space in each site, InD. This clearly puts 
all random singlet phases on the same footing even when 
arising in different microscopic systems. In addition, one 
may ask whether such a redefined measure may always 
reduce along RG flows, even those connecting pure fixed 
points to random ones. 

An interesting difference between the spin-1 case we 
study and previous work, is that for RS phases the ex- 
act logarithmic divergence of entanglement entropy can 
be found in closed form using the RG history approach, 
while for spin-1, the true entanglement entropy seems to 
depend on the entanglement values of an infinite num- 
ber of nonequivalent, irreducible sub-chains. This turn 
of events motivated us to define the 'reduced entropy' 
in Sec. V. This simplified measure of the entanglement 
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counts directly how many SBS's connect the segment L 
with the rest of the chain, and neglects the correlation 
induced by the SBS's. Our analysis allowed an exact 
calculation of the reduced entanglement, which yielded: 



2S 

^reduced ^ 1 3 « 1.329. 

■' 12 



(151) 



A numerical test using a computerized application of the 
RSRG confirmed this exact analj^ical result (see Sec. 
VF), thus affirming the history-dependence segment of 
our approach. 

In light of the simplicity of the definition of the reduced 
entanglement, it is interesting to ask how is relates to the 
real quantum entanglement. Since it neglects correlations 
between singlets, the reduced entanglement is most likely 
an overestimate of the quantum entanglement. This re- 
mains to be proved generally. In the spin- 1/2 case, the 
reduced entanglement and the real entanglement coin- 
cide. For the spin-1 RS-Haldane critical point, we have: 



^reduced 
^eff 



1.329 > 1.232 w c 



Ire) 
■eff 



(152) 



Hence the reduced entropy is bigger by about 10%. If 
there were some simplifying relation in the entanglement 
of these irreducible units, as there is for the reduced en- 
tanglement, then our method would give a closed form, 
but barring that it seems that entanglement entropy is 
a less natural object in the real-space renormalization 
group for higher spin than the (unphysical) reduced en- 
tropy. 

Future directions include, needless to say, the calcu- 
lation of the universal entanglement entropy in other 
infinite-randomness non RS fixed points. Our method 
can be straightforwardly applied to s > 1 chains. 
But an exciting possibility is the consideration of the 
infinite-randomness fixed points arising in non-abelian 
spin chains. The random-singlet subset of this class of 
fixed points were recently studied in Ref. 19. 
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APPENDIX A: CALCULATION OF Nu, N- 



n 



In Sec. VI B we defined the states: 



IT (1, 1) T) = 44 "n (aVbl+i - bl&l+i) |0) 

IT (1,1) i) = 4&rn - io) 

li (1, 1) T) = h\al "n - |0) 

li (1,1) i) = stst"n (aRV-^l«l+i)io)- 



(Al) 



i=0 



We found that these states do not constitute a orthogonal 
set since: 



(T (1,1) ill (1,1) T) = (-i)". 



(A2) 



To be able to carry out the procedure outlined in Sec. 

VI B we need to calculate the norms of these states, which 
will than allow a transformation into an orthonormal 
combination. 

From standard transfer matrix techniques, and Wick 
contractions, we find for N-^y. 
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